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We studied the evolution of isolated strange stars (SSs) synthetically, considering 
the influence of r-mode instability. Our results show that the cooling of SSs with non- 
ultrastrong magnetic fields is delayed by heating due to r-mode damping for millions of 
years, while the spin-down of the stars is dominated by gravitational radiation (GR). 
Especially for the SSs in a possible existing colour-flavour locked (CFL) phase, the 
effect of r-mode instability on the evolution of stars becomes extremely important 
because the viscosity, neutrino emissivity and specific heat involving pairing quarks 
are blocked. It leads to the cooling of these colour superconducting stars being very 
slow and the stars can remain at high temperature for millions of years, which differs 
completely from previous understanding. In this case, an SS in CFL phase can be 
located at the bottom of its r-mode instability window for a long time, but does not 
spin-down to a very low frequency for hours. 
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1 INTRODUCTION 

The r-modes in a perfect fluid star succumb to gravitational 
radiation (GR) driven Chandrasekhar-Friedman-Schutz in- 
stability for all rates of stellar rotation, and they arise due 
to the action of the Coriolis force with positive feedback in- 
creasing GR (Andersson 1998; Friedman & Morsink 1998). 
However, actually, the viscosity of stellar matter can hold 
back the growth of the modes effectively. Based on the 
competition between the destabilizing effect of GR and the 
damping effect of viscosity, Owen et al. (1998) modelled the 
evolution of the r-modes and described the spin-down of neu- 
tron stars as losing angular momentum to GR. Afterwards, 
Ho & Lai (2000) improved the model considering the conser- 
vation of wave action and containing the magnetic braking 
due to magnetic dipole radiation (MDR). 

The viscosity of stellar matter is decisive for the r-mode 
evolution and hence for the spin-down and cooling of the 
star. The bulk viscosity of strange quark matter in the nor- 
mal phase is several orders of magnitude stronger than that 
of normal nuclear matter. As a consequence, the instabil- 
ity in a strange star (SS) can only be driven at relatively 
lower temperature (< 10** K), which means not very early 
ages, differing from the case of a neutron star in which the 
instability acts immediately at the birth of the star (Owen 
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et al. 1998; Watts & Andersson 2002). Thus Madsen (2000) 
argued the evolution of SSs can be very different from that 
of neutron stars. Andersson, Jones & Kokkotas (2002) cal- 
culated the evolution of SSs thoroughly using the model of 
Ho & Lai (2000) (but without a magnetic field) in the pres- 
ence of accretion to point out that SSs may be persistent 
sources of GR. However, most works always focused on the 
stars with r-mode instability as sources of GR, although the 
authors also claimed that the r-modes could change the spin 
evolution of the stars due to GR and induce heating due to 
viscous dissipation of the modes. 

We also note the work of Gusakov, Yakovlev, & Gnedin 
(2005) who studied the thermal evolution of a pulsating neu- 
tron star, taking into account the pulsation damping due to 
viscosity and accompanying heating of the star. This idea 
leads us to pay attention to the effect of the r-mode insta- 
bility on the evolution of compact stars. As shown in fig. 
5 in Watts & Andersson (2002), by heating due to r-mode 
damping, the cooling of a nascent neutron star is delayed in 
the first few decades. In addition, the calculation in Ander- 
sson et al. (2002) implies that the history of the r-modes in 
isolated SSs could be longer than neutron stars. According 
to the above works, we consider that a careful investigation 
of the effect of the r-mode instability on the evolution of SSs 
is an interesting topic, which we will focus on in this paper. 
As expected, our results show that, although the r-modes 
may be too small to induce a detectable GR at the age of a 
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million years, the accompanying heating effect still can not 
be ignored in SSs. 

In particular, our attention is attracted by the same 
consideration of colour superconducting SSs. For strange 
quark matter at a sufficiently high density, phenomenologi- 
cal and microscopic studies predict it could undergo a phase 
transition into a colour superconducting state, such as the 
typical cases of two-ffavour colour superconductivity (2SC) 
and the colour-flavour locked (CFL) phase (Shovkovy 2005; 
Afford 2004). Theoretical approaches also concur that the 
pairing gap (A) in the quark spectrum is approximately 
~ lOOMeV for baryon densities existing in the interior of 
compact stars. In this case, the reaction rates involving two 
quarks are reduced by a factor exp(— 2A/fcBr), assuming 
equal behaviour for all quark flavours. Because the quarks 
in the 2SC phase pair in part, the corresponding SSs will 
have similar behaviours to SSs in normal phase (NSS here- 
after; Blaschke, Klahn, & Voskresensky 2000; Madsen 2000; 
Yu & Zheng 2006). Therefore, here we only account for SSs 
in CFL phase (CSS hereafter) with A = lOOMeV, in which 
all flavour quarks pair. As a consequence, all of the viscos- 
ity, the neutrino emissivity and the quark speciflc heat of 
CSSs are reduced. Thus the effect of the r-mode instabil- 
ity on the evolution of CSSs may become more important. 
As shown in Section 5, the dramatic variation in the stellar 
evolution changes our previous understanding of CSSs that 
the stars are very cold and slow revolving objects (Blaschke 
et al 2000; Madsen 2000). 

This paper is organized as follows: In Section 2 and 3 we 
recall the r-mode instability and the model of spin evolution 
due to gravitational and magnetic braking respectively. The 
thermal evolution equation of SSs considering the heating 
due to the viscous dissipation of the r-modes is given in 
Section 4. Section 5 presents the numerical results of the 
evolution of NSSs and CSSs. Finally, we give the conclusion 
and our discussions in the last section. 



represent large-scale oscillating currents that move (approxi- 
mately) along the equipotential surfaces of the rotating star. 
The density perturbation associated with the r-modes can 
be deduced by evaluating the inner product of v (the unper- 
turbed fluid velocity) with the perturbed Euler equation, 
and the equation for the perturbed gravitational potential 
(Lindblom et ale 1998): 
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The quantity S'i' is proportional to the perturbed gravita- 
tional potential 5$ and is the solution to the ordinary dif- 
ferential equation 



d^S-^(r) 2 d(5^(r) 

rf75 r - 



+ [4-Gpf 




-1)('+2) 1 
"T2 



5*(r) 



dp ;' r \'+l 
l + lPdp 



i^y 



(5) 



The r-modes evolve with time dependence e'^'^^'' ^" as 
a consequence of ordinary hydrodynamics and the influence 
of the various dissipative processes. The real part of the fre- 
quency of these modes, lj, is given in equation (3), while 
the imaginary part 1/r is determined by the effects of GR, 
viscosity, etc. The simplest way to evaluate 1/r is to com- 
pute the time derivative of the energy E of the mode (as 
measured in the rotating frame). E can be expressed as a 
real quadratic functional of the fluid perturbations: 
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Thus, the time derivative of E is related to the imaginary 
part of the frequency 1/r, which can be conveniently de- 
composed into two parts associated with GR and viscosity, 
respectively. It reads 
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2 THE B-MODE INSTABILITY 

The r-modes of rotating Newtonian stars are generally de- 
flned to be solutions of the perturbed fluid equations hav- 
ing (Eulerian) velocity perturbations of the form (Lindblom, 
Owen & Morsink 1998) 



Sv = aRQ, 
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where R and il are the radius and angular velocity of the 
unperturbed star, a is an arbitrary constant considered as 
the amplitude of the r-modes (Owen et ale 1998), and F,^ 
is the magnetic-type vector spherical harmonic deflned by 

Yi'L = [l{l + 1)]"'/VV X (rVYim). (2) 

Papaloizou & Pringle (1978) flrst showed that the Euler 
equation for r-modes determines the frequencies as 
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Further use of the Euler equation (as flrst noted by Provost, 
Berthomieu & Rocca 1981) in the barotropic case deter- 
mines that only I — m r-modes exist and that 5v must have 
the radial dependence given by equation (1). These modes 



For n = 1 polytrope, the time-scale for GR is calculated by 
Lindblom, MendeU & Owen (1999) 
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where p is the mean density of the star. The viscous time- 
scale is contributed by shear and bulk viscosity as 



T"v = (rsv + Tb 



(9) 



For strange quark matter, we can give (Lindblom et al. 1999; 
Madsen 2000) 



r,v = 5.41 X W^'al^l-.T^/^ 
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where Qc.o.i, mioo, and Tg are the strong coupling in units 
of 0.1, the mass of strange quark in units of lOOMeV and 
the interior temperature of the stars in units of lO^K, re- 
spectively. When these viscosities due to quark reactions 
are blocked in the CFL phase, the dissipation of the r-modes 
will be mainly dominated by a smaller shear viscosity due to 
electron-electron scattering. The corresponding time-scale is 
(Madsen 2000) 
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From equation (7), it can be seen that the r-modes are un- 
stable if Tg"^+T\r^ < 0. In this case, a small perturbation will 
lead to an unbounded growth of the modes. The competition 
between the destabilizing effect of GR that is dependent on 
spin frequency and the damping effect of the temperature- 
dependent viscosity gives an instability window in the v — T 
plane, as shown by the shadows in the following figures. We 
can see that, compared with the window of an NSS, the 
window of a CSS expands significantly because its viscosity 
involving quarks are blocked. 



3 SPIN EVOLUTION 

We employ a simple phenomenological spin evolution model 
proposed by Ho & Lai (2000), which is analogous to that 
devised by Owen et al. (1998). In this case, the total angular 
momentum of a star can be obtained as the sum of the bulk 
angular momentum and the canonical angular momentum 
of the r-modes Jc 



J^m + Je, 



(13) 



where / = IMR^ [I = 0.261 for n — 1 polytrope) represents 
the moment of inertia of the star. The canonical angular 
momentum of an r-mode can be expressed in terms of the 
velocity perturbation 5v by (Friedman & Schutz 1978) 
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For the (dominant) I = m = 2 current multipole of the r- 
modes, the above expression reduces (at the lowest order in 
n) to 



Jc = -'^a'^JMR^n, 
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where J — 1.635 x 10^^ for n — 1 polytrope. The canonical 
angular momentum of the mode can increase through GR 
back reaction and decrease by transferring angular momen- 
tum to the star through viscosity: 

(Jbo c ^'J C ^'J C 

dt Tu Tv 
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On the other hand, adding the reverse angular momentum 
transferred from the r-mode, the bulk angular momentum of 
the star decreases. At the same time, the star also undergoes 
braking by MDR. Thus the rotation of the star is determined 
by 
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where Tm = 1.69 x 10^ B-^^2i^/V''^Gp) ^s is the magnetic 
braking time-scale, and B12 is the magnetic field intensity 
in units of 10^^ G. 

Submitting equation (15) into equations (16) and (17), 
we can obtain the coupled evolution equations of the ampli- 
tude a of the r-mode and the angular velocity fl: 
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where Q = 3J/2/ — 0.094. In some cases, the small initial 
perturbation could increase to a large value so that non- 
linear effects can no longer be ignored. As postulated by 
Owen et al. (1998), there may exist a saturation amplitude 
Q^ — K. When a^ > k, the growth of the amplitude of the 
mode stops, then 
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In this stage, the spin evolution will be described as 
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Of fundamental importance in judging the astrophysical rel- 
evance of the instability is the determination of the satura- 
tion amplitude of the mode. Lindblom, Tohline & Vallisneri 
(2001, 2002) suggested that the non-linear saturation ampli- 
tude may be set by dissipation of energy in the production of 
shock waves. However, Gressman et al. (2002) argued that 
the decay of the amplitude of the order unity is due to leak- 
ing of energy into other fluid modes, leading to a differential 
rotation configuration. Afterward, the coupling between the 
r-modes and other modes was analyzed (Arras et ale 2003). 
On the other hand, the role of differential rotation in the 
evolution of the r-mode was also studied thoroughly (Sa & 
Tome 2005). All these works obtained a credible conclusion 
that the saturation amplitude of the r-mode may be not 
larger than the small value of 10~^. Thus, we take k — 10"^ 
in our calculations. The saturation amplitude is important 
for the evolution of CSSs because the r-mode in CSSs will 
retain the saturated state for a long time. However, to a cer- 
tain extent, the uncertainty of the saturation value can only 
influence our numerical results quantitatively, but cannot 
change the conclusions essentially. Note that, for NSSs, this 
artificial cut of the amplitude is needless, as is also claimed 
by Andersson et al. (2002). 



4 THERMAL EVOLUTION 

An SS can be divided roughly into an inner, approximately 
isothermal, quark core and a thin, solid, nuclear crust (Al- 
cock, Farhi & Olinto 1986). For the relationship between the 
interior temperature T and the surface temperature Ts, we 
apply the result given by Gudmundsson, Pethick & Epstein 
(1983), which is valid for a crust with the density at the base 
that is just larger than 10*g cm^'\ Thus, the emissivity due 
to the surface photon emission can be written as 

S.hoton = 1.24xl0''r9''' erg cm"' s"\ 



^photon 



(22) 



Besides the contribution from this thermal radiation, SSs 
also cool via internal neutrino emission. In strange quark 
matter, there are three main neutrino processes: (i) di- 
rect Urea processes d^uei' and ue^>dv, (ii) modified 
Urea processes dq^mqeu and uqe~^dqp; (iii) the quark 
bremsstrahlung processes qiq2^qiq2i^i'- The corresponding 
neutrino emissivities read respectively (Iwamoto, 1982) 
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e'^'^\ . = 2.1 

neutrino 
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where pb is the baryon number density and po ~ 0.17 fm~ 
is the nuclear saturation density. Ye = pe/pb is the elec- 
tron fraction. On the other hand, SSs also can be heated. 
During GR, in addition to the losing of angular momen- 
tum, the r-modes also can lose energy via GR and neu- 
trino emission (from bulk viscosity) and can deposit energy 
into the thermal state of the star via shear viscous dissipa- 
tion (Owen et al. 1998). As defined in equation (6), for the 
Z = m = 2 current multipole, the energy of the r-mode is 
E = \a^JMR^Q?, and the heating rate is 
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Thus the thermal evolution equation is written as 
dT 

^ , — -^neutrino -^^photon i -ngv, V j 

where -Lnoutrino and Lphoton are the neutrino luminosity and 
the surface photon luminosity respectively, and C is the total 
specific heat, which is mainly contributed by quarks and 
electrons (Blaschke et al. 2000): 



Ca~2.5xl0^ 



Co~0.6xl0^ 



2/3 



Tg erg cm K 



Vcpb 



2/3 



Tg erg cm K 



(28) 



(29) 



In the calculation for CSSs, the reductive factor due to the 
pairing gap should be added in equations (23)-(25) and (28). 
In other words, the electron specific heat, the surface photon 
emission and the heating term become more important for 
the thermal evolution of CSSs. 
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Figure 1. The evolution of NSSs (with canonical parameters 
M = IAMq and i? = 10 km) for different magnetic fields. The 
labels for the open cycles represent the value of log(t yr^-'^). 
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5 EVOLUTION CURVES 

In our calculations, wo take the initial temperature To = 
lO^^K, the initial angular velocity fio = S^yrGp, and a rep- 
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resentative set of parameters of Y^ 
lOOMeV for SSs.^ 

We firstly discuss the evolution of NSSs. The stellar evo- 
lution tracks in the f — T plane for different magnetic fields 
are shown in Fig. 1 and several typical ages are marked on 
the tracks using open cycles. Fig. 2 shows the amplitude evo- 
lution of the r-mode with time corresponding to the same 
field cases. We can see from Fig. 1(b), (c) and (d), just as 
calculated by Andersson et al. (2002) for B = Q, that the 
stars for non-ultrastrong {B < 10^'"* G) field cases will reach 

^ We provide a brief review of the relevant properties of strange 
quark matter in the previous sections in order to discuss the evo- 
lution characteristic of SSs in this section. However, a more com- 
prehensive investigation of these properties is needed to compare 
the theoretical curve with the observational data. Comments on 
an early version of this paper helped us to note that the con- 
tribution to viscosity in the CFL phase from superfluid phonons 
may be more important than that of electrons (Manuel, Dobada 
& Llanes-Estrada 2005), although the physical inputs left un- 
mentioned in this paper may not affect the general astronomical 
picture. 
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Figure 2. The evolution of the amplitude of the r-mode in NSSs 
for different magnetic fields. 



the instability windows several days after birth as a result of 
temperature decrease and will then evolve along the bound- 
ary of the windows. However, we also find that the specific 
history of the stars is sensitively dependent on the magnetic 
field intensity. It can be seen from the labelled ages that, the 
weaker the field, the slower the spin-down and the cooling. 
Especially with very weak {B < lO^'^G) fields, the stars can 
remain at a very high spin frequency, such as 400Hz, even at 
the age of over a million years. Finally, the evolution tracks 
will depart from the instability windows at the moment that 
the magnetic braking effect exceeds that of GR (see the up- 
per panel of Fig. 3). Thus, the stronger the field, the earlier 
the departure. Synchronously, this departure also leads to 
the end (a catastrophically decay) of the r-mode as shown 
in Fig. 2. Before this abrupt drop, however, the variation of 
the mode is small apart from the first several decades, during 
which the amplitude of the r-mode oscillates greatly until it 
converges to a value of ~ 10^^ as discussed by Andersson 
et al. (2002). At the extreme, a star with an ultrastrong 
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Figure 3. The upper panel shows the spin evolution of NSSs 
for different magnetic fields: B = 10^"' G (dash-dotted curve); 
B = IQl^G (dotted curve); B = 10"G (sohd curve); B = IQIOG 
(dashed curve); and the thin solid curves represent the spin-dovifn 
induced by MDR purely corresponding to the same field cases 
from left to right. The lower panel shows the surface tempera- 
ture evolution of NSSs for the cases in the upper panel, and the 
thin solid curve corresponds to cooling in the absence of r-mode 
instability. 
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Figure 4. The evolution of CSSs (with canonical parameters 
M = IAMq and i? = 10 km) for different magnetic fields. The 
labels for the open cycles represent the value of log(i yr~^). 
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Figure 5. The evolution of the amplitude of the r-mode in CSSs 
for different magnetic fields. 




Figure 6. The upper panel shows the spin evolution of CSSs 
for different magnetic fields: B = 10^'^G (dash-dotted curve); 
B = IQl^G (dotted curve); B = 10"G (sohd curve); B = IO^OG 
(dashed curve); and the thin solid curves represent the spin-down 
induced by MDR purely corresponding to the same field cases 
from left to right. The lower panel shows the surface tempera- 
ture evolution of CSSs for the cases in the upper panel, and the 
thin solid curve corresponds to cooling in the absence of r-mode 
instability. 
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{B > 10^^ G) field decelerates so rapidly due to MDR that 
the track will not come into contact with the instability win- 
dow as shown in Fig. 1(a). Thus, no r-mode arises in the star 
throughout its life as shown in Fig. 2(a). To be clear, Fig. 3 
shows the evolution curves of spin frequency (upper panel) 
and surface temperature (lower panel) corresponding to the 
cases shown in Fig. 1. Not surprising, being governed by 
Eq. (19), the spin history of the stars with non-ultrastrong 
(B < 10^*^0) fields can be divided into two stages: (i) the 
stage due to angular momentum transfer by viscosity (los- 
ing angular momentum via GR ultimately) at early ages in 
the presence of the r-mode; and (ii) the stage due to MDR 
in the absence of the r-mode. On the other hand, as shown 
in the lower panel, the cooling of the stars with weak fields 
is delayed under the impact of heating due to the viscous 
dissipation of the r-mode, whereas the curve of the strong 
field superposes the curve without the consideration of the 
r-mode instability. 

In the second situation, colour superconductivity sup- 
presses the cooling effect extremely, whereas the heating ef- 
fect is enhanced considerably. The evolution tracks of CSSs 
in the v — T plane for different magnetic fields and the cor- 
responding amplitude evolution curves of the r-mode are 
shown in Fig. 4 and 5, respectively. The same as Figs 3 but 
for CSSs are also shown in Fig. 6. Differing from NSSs, the 
r-mode in CSSs can appear at the birth of the stars and 
quickly (~ 1000s, as shown clearly in Fig. 5) rise to the sat- 
urated state because the bulk viscosity is blocked. However, 
at the same time, the temperature nearly does not decrease 
because the neutrino cooling effect is suppressed. However, 
because the specific heat is only contributed from electrons, 
the temperature decreases finally in the following several 
hundred years until the photon luminosity is nearly com- 
pensated by the heat gain from the shear viscous dissipation 
of the r-mode. From then on, the stars should experience a 
long-term slow cooling stage. As discussed by Blaschke et al. 
(2000) and Yu & Zheng (2006) previously, CSSs will become 
very cold at early ages (< lOOOyrs) if the heating effect is 
ignored. However, the presence of the long-term slow cooling 
stage can change the situation completely as shown in Fig. 
6 (lower panel). We can see the surface temperature can re- 
main at a high value of ~ 10® K lasting several million years. 
The spin-down of CSSs that is shown in the upper panel of 
Fig. 6 is similar to the one of NSSs, but the stage dominated 
by GR is much longer because the life of the r-mode in CSSs 
is longer as shown in Fig. 5. 



6 CONCLUSION AND DISCUSSIONS 

We have studied the evolution of SSs synthetically, consid- 
ering the influence of r-mode instability. The spin history is 
divided into two stages and the heating effect due to viscous 
dissipation of the r-mode can delay the cooling of SSs with 
non-ultrastrong fields remarkably. Our prime objective is to 
investigate the behaviour of SSs in colour superconducting 
phase (to be specific, the CFL phase). This problem is es- 
timated simply by Madsen (2000), who argued that CSSs 
should spin-down to a very low spin frequency for hours due 
to the r-mode instability. Conversely, our results indicate 
that the stars can be located at the bottom of the window 
for millions of years due to their characteristic cooling be- 



haviour. A CSS could be a rapid revolving pulsar. Differing 
from NSSs, the cooling of CSSs is changed completely by 
heating due to r-mode damping. Except for the first several 
hundred years, CSSs cool very slowly. In other words, the 
stars can remain at a high temperature over several million 
years but do not become very cold at early ages as calculated 
in the absence of the heating effect. Thus, we believe that 
it is possible to detect a hot source being a CSS but not, as 
argued by Blaschke et al. (2000), that CSSs are rejected by 
the X-ray data. 

In our opinion, when we talk about stellar evolution, we 
must consider the cooling and the spin-down synthetically 
because the rotation of the star is a large possible energy 
store for the heating of the star. We note that Reisenegger 
(1995) proposed a chemical heating mechanism due to the 
derivative from the /3— equilibrium of neutron matter with 
the spin-down of neutron stars. The same issue was dis- 
cussed in NSSs by Cheng & Dai (1996), who presented that 
chemical heating delays the later (> lO^yrs) cooling of NSSs 
significantly. In this paper, however, it is the r-mode insta- 
bility that couples the cooling and the spin-down: the kinetic 
energy transforms into the canonical energy of the r-modes 
via GR, then is deposited into a thermal state via the dissi- 
pation of the r-modes. Our results show that the heating due 
to the r-mode damping delays the earlier (< lO^yrs) cool- 
ing remarkably. To summarize, we argue that the heating 
due to spin-down plays an important role in the evolution 
of an NSS throughout its life. On the other hand, for a CSS, 
the chemical heating is blocked because the reactions be- 
tween the quarks are suppressed by the pairing gap. Thus, 
the heating due to the r-mode damping becomes more de- 
cisive for its evolution. As discussed above, the spin-down 
can impact the cooling of an SS as an energy source. Besides 
this direct influence, the spin evolution also can determine 
the change of other properties of the star, such as the vari- 
ation of the density of the star. We therefore think that the 
spin-down can also impact the cooling behaviour indirectly. 
Looking at more detailed calculations, previous modelling of 
thermal evolution of SSs based on the assumption that the 
spin-down is only determined by MDR needs to be recon- 
sidered because it is GR but not MDR that dominates the 
flrst spin-down of the star in the presence of the r-modes. 

Using NSSs, several authors (Madsen, 2000; Andorsson 
et al. 2002) try to explain the clustering of low mass X-ray 
binaries in the v — T plane. Our results (Fig. 1) show that, 
for weak-field {B < 10^" G) cases, MDR has no influence on 
the spin-down of the stars that are younger than 10 million 
years old. This implies the existence of a weak magnetic field 
in millisecond pulsars if we follow the same philosophy that 
has been used by Madsen (2000). However, in our opinion, 
the final clarification of this problem should dependent on 
careful consideration of the accretion of stars as Watts & 
Andersson (2002) did for neutron stars. 
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